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We calculate the dimensionless Fermi liquid parameters (FLPs), Fq^" 1 and F^ 171 , 
for spin asymmetric dense quark matter. In general, the FLPs are infrared divergent 
due to the exchange of massless gluons. To remove such divergences, the Hard 
Density Loop (HDL) corrected gluon propagator is used. The FLPs so determined 
are then invoked to calculate magnetic properties such as magnetization (M) and 
magnetic susceptibility xm of spin polarized quark matter. Finally, we investigate 
the possibility of magnetic instability by studying the density dependence of (M) 
and xm- 
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I. INTRODUCTION 



The study of strongly interacting matter has been an area of contemporary research for 
quite sometime now. Such studies are usually made in the extreme condition of temperature 
and/or density. The high temperature (T) studies are more relevant to the ultra-relativistic 
heavy ion collisions while the investigations involving high chemical potential (/i) or extreme 
case of cold matter are more appropriate to astrophysics It should, however, be noted 
that efforts are being directed recently to also study the properties of a very dense system in 
;he laboratory where matter with predominantly large chemical potential might be formed 
2|. We here restrict ourselves to zero temperature and investigate some of the properties of 
quark matter in presence of a weak magnetic field. 

It has been shown recently that the degenerate quark matter can show para-ferro phase 
transition below a critical density {3]. To examine this possibility, in Ref. |3|, a variational 
calculation was performed. Subsequently, various other calculations were also performed in 



different formalisms to investigate such a possibility with varied conclusions [4HL1]. 

The issue of spontaneous phase transition in dense quark system at zero temperature 
was also examined in 12j by invoking Relativistic Fermi Liquid Theory (RFLT). In partic- 



ular, this was accomplished by calculating the chemical potential (//) and energy density of 
degenerate quark matter in terms of the Landau Parameters (LPs). The RFLT was first 
developed by Baym and Chin [13I . to study the properties of high density nuclear mat- 



ter. However, the formalism developed in Ref.[13| is valid for unpolarized matter and LPs 
calculated there are spin averaged. Here, on the other hand, we deal with polarized quark 
matter which requires evaluation of the LPs with explicit spin dependencies. 



Recently, in [15l . Il6j the authors have studied the magnetic properties of degenerate 
quark matter in presence of weak uniform external magnetic field B. Similar investigation 
was also made in Ref. 11] by evaluating the effective potential and employing quark magnetic 
moment as an order parameter. These calculations were, however, restricted to the case of 
unpolarized matter. On the contrary, our concern here is the magnetic properties of polarized 
quark system. Consequently, we first determine various spin combination of LPs such as 
spin symmetric (F^i ' sym ) and spin anti-symmetric (F+(~*)> asyms ) parameters and express 
quantities like magnetization and magnetic susceptibility in terms of these parameters. It is 



needless to mention that unlike 



11 



15 



161 ] . the expressions for xm and (M), as presented 
here, depend on the spin polarization parameter £ = (n+ — n~)/(n+ +n~), where n+ and 
n~ correspond to densities of spin-up and spin-down quarks, respectively. 

It is well known that the calculations of LPs require evaluation of the forward scatter- 
ing amplitudes which are plagued with infrared divergences arising out of the exchange of 
massless gluons. Formally, such divergences can be removed by using HDL corrected gluon 
propagator. This can also be achieved by introducing screening mass for the gluons. Such 
regularizations are necessary for the evaluation of individual LPs. On the other hand, in 
various physical quantities like the ones we calculate here, the LPs appear in particular com- 
binations where such divergences cancel at least to the order with which we are presently 
concerned. 

The plan of the article is as follows. In Sec. II we derive the expressions of LPs for 
polarized quark matter. In Sec III, we calculate magnetic susceptibility in terms of LPs 
with explicit spin dependencies both with bare and HDL corrected gluon propagator. In 
Sec. IV we summarize and conclude. 



3 



II. SYMMETRIC AND ANTI-SYMMETRIC LANDAU PARAMETERS 

In this section we calculate LPs for spin polarized quark matter. We are dealing with 
quasi-particles whose spins are all eigenstates of the spin along a given direction viz. z. The 
quasiparticle interaction can be written as the sum of two parts viz. spin symmetric {fppT 1 ) 
and anti-symmetric {fppf m ) parameters [l^, 16| : 



f ss' _ fsym , / A rasym 
Jpp' ~ Jpp' v ' J Jpp' 



(1) 



Assuming that the spins are randomly oriented with respect to the momentum, we take 
average over the angles Q\ and 82 corresponding to spins s and s'. The angular averaged 
interaction parameter is given by 12[: 



fss> 

J pp 1 



p=Pf,p-p f 



2 . f ss> 

J mi' 



P=P f ,P-Pf 



(2) 



An Jpp 

1 Here the spin may be either parallel (s = s') or anti-parallel (s = —s') jsl, [l^. Thus the 
scattering possibilities are denoted by (+, +), (+, — ), (— , — ) etc. The interaction parameters 
can now be redefined as, 



f++ _ rsym , rasym _ r 

J pp' J pp' ' J pp' J pp' 



J pp' 



psym rasym e — y 

J pp' J pp' J pp' 



(3) 



Once these interaction parameters are known, the FLPs can be determined by expanding 
fpp, into the Legendre polynomial: 



r' = (2/ + i) J ^Pi(cos6)f;;'„ 



(4) 



where cosO = p ■ p' . We define symmetric and anti-symmetr ic p art of LPs y*'^" 1 ^^" 1 ) w hat 
one does to dealing with the isospins in nuclear matter 12|, |l4| : 



[1] denoted hereafter as /if,' = /if, 



4 



fi 



fi 



+(-),sym 
I 

+{—),asyra 



++(~) 



+-(-+) 



+-(-+) 



(5) 



It should be noted here that, f^ p , = f vx }~ . 



VP 



The dimensionless LPs are defined as p^ s y m ( as y m ) = jys^js,syrn{asym) w Y iere j\r s (0) 
is the density of states at the Fermi surface, which can be written as, 



N s (0) 



d 3 p s , 
^d(e ps - /I J 



(27T) 

N cP f ( dp 



2tt 2 



de 



(6) 



Here, iV c is the color factor, e ps and /i s are the spin dependent quasi-particle energy and 
chemical potential respectively. It is evident from Eq.flE]) that for spin polarized matter, the 
density of states is spin dependent. This, as we shall see, makes the calculation cumbersome. 



In the above expression (dp/de 



by [12|,|l5l 



ps , 



p=pf 



is the inverse Fermi velocity (1/vj), where v s f is given 
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p) N c p s f fi' g y m 



2tt 2 



(7) 



With the bare propagator, the angular averaged spin dependent interaction parameter 
yields [12] 



Jp P > 



p=p'=pJ 



9sj 2 pj 2 (l-cosl 



2mi 



COS I 



2m q p^ 2 
3(£j + m q 



f + 7 
J pp' 



p=p} ,p =p f 



+2 

m q pj 



in 



+ 



3(e~j + m q ) 3(e f + m q ) 



1 



£f£f + PfPf cosl 



(9) 



Here, m q is the quark mass, p^ = p/(l ± £) 1//3 , = (p^ 2 + m 2 ) 1 / 2 and pj is the Fermi 
momentum of the unpolarized matter (£ = 0). Similarly, f~7 can be obtained by replacing 
with pj and with ej in Eq.flH]). One can find dimensionless LPs, F^ 71 and jP^ m 
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(suppressing spin indices) by considering OGE interaction. But both of these (_p^ m ( a ^ m )) 
exhibit infrared divergences because of the term (1 — cos 9) that appear in the denominator of 
the interaction parameter (see Eq. (15|) ). This divergence d isap pears if one uses HDL corrected 
gluon propagator to evaluate the scattering amplitudes [IT] - 

To construct HDL corrected gluon propagator with explicit spin dependence one needs to 
evaluate the expressions for longitudinal (n^) and transverse (Ht) polarization which have 
been derived in We borrow the results directly: 



n L (k ,k) 

n T (k ,k) 



3_ 

Air 2 

g 2 



(cg-i)E^ 



s=± 




16n 



s=± 



2C n 



1 



In 



Co 



f 



C - v' f 



(10) 
(11) 



"f 

Here, Co = fco/|A;|, is the dimensionless variable and = p^/e^. It might be noted here, 
that the expressions for and look rather similar to what one obtains in the case of 
unpolarized matter (£ = 0) 



with only difference in v f . In the static limit i.e. Co 0, 



the spin dependent Debye mass (mo) is given by 



n. 



m 2 D 



JL 
An 2 



(12) 



s=± 



It is to be mentioned here, that the screening mass of the gluon is spin dependent and 



the transverse gluons are screened only dynamically 
combination of dimensionless LPs are found to be 



15 



161 ] . With these, the symmetric 



rp+,sym 
r 



g Pf r l 
144vr 2 I e 



+- 



/ 

12 + 



12 - 



\2m\ + 12m 2 q ej + 3m q m 2 D + Apj 2 m q + 3m 2 D ej / Apf 



P~f 2 (m q + sy 



PfPf 



ej)(m q + e f 



-m q [3m 



2 (f + 

D\ £ f 



2(e+pj 2 + ejpf)] + 3m 2 D e+ej} 



x In 



/2m 2 -m 2 D + Zp^fPf - ^ £+ f e ) 
\2m 2 — m 2 D — 2p~jpJ — 2e~je~ } 



In 1 + 



m 



D 



m 2 q [3m 2 D -2{pf+p- 2 )} 



(13) 
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,sym 



9 2 p1 r 1 



X 



X 



48vr 2 I e 



I 



12m q + 12m 2 e^ + 3m q m 2 D + 4p~j: 2 m q + 3m 2 D e^ 
p+ 2 {m q + e+) 



m 



D 



ln ( 1 + J:) 



2 P/ ' '"£) 



5 



m3[3rr& - 2(p+ 2 + pj 2 )} + mjSm^e} + ej) - 2{e+pf + ejpj -2 )] + 3m 2 D e+e) 



PfPf ijn q + e^){m q + e f ) 



f 

2m 2 q -m 2 D -2eje] ^ ^ ( 2m l q - m 2 D + 2pjpJ - 2ep' f 



2p fPf 



2m 2 — m D — 2pjPj- — 2e^ e ^ 



(14) 
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FIG. 1: Dimensionless LPs as a function of Fermi momentum for unpolarized and polarized quark 
matter. Symmetric and anti-symmetric combination of LPs are plotted in (a) and (b) respectively. 



In deriving Eqs. (ll3[) and (III"]) , we consider exchange of longitudinal gluons only. In 
Eqs. (I13j) and (1141) . the term in the first square bracket arises due to the scattering of like- 
spin states (++), while the latter comes from the scattering of unlike-spin states (-| — ). 
Similarly one may determine other combination of LPs like FQ{ sym , p+^ as y m ^ p-^ s v m e ^- c _ j n 
Fig{H density dependence of symmetric and anti-symmetric combination of dimensionless 
LPs is shown. Similar plots for the LPs in isospin asymmetric nuclear matter can be found 



m [19]. There, however, the calculated LPs are finite, as the nucleon-nucleon interactions 
involve exchanges of massive mesons like a, u, 5 and p etc. It is interesting to note that the 
results of isospin asymmetric nuclear matter for the LPs are qualitatively same as those of 
dense quark system. 
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III. MAGNETIC SUSCEPTIBILITY 



Now, we proceed to calculate the magnetic susceptibility for which an uni 
field B is applied along the z axis. The magnetic susceptibility is defined as 



orm magnetic 
15, 



16] 



Xm 



E 



d(M) j 



dB 



B=0 



(15) 



where (M) f is the magnetization for each flavor. Here, the magnetic field is considered to 
be significantly weak for which the spinors remain unaffected and only modification enters 
through the single particle energy. Here, we consider one flavor quark matter and suppress 
the flavor indices. 

In presence of constant magnetic field B, the magnetization depends on the difference of 
the number densities 5np Sym = 5n p , s= i — 5n p ^ s= ^i, where 



n ps = [1 + exp (3(e ps - /i - -g D {p)fi q sB)] 1 . 



(16) 



In the last equation, fi q denotes the Dirac magnetic moment and go{p) is the gyromagnetic 
ratio. The magnetization is given by [jjj] 



(M) = T N °J WT Mp)5n7m 



(17) 



For constant magnetic field, the variation of the distribution function yields 15|, 1 161 ] . 



5n 



5n 



ps 



ps 



5e 



ps 



1 f ; 

-g D (p)v q sB + N C J2 J j^yf s P p> 6n P's' 



and <5rip SJ/m is therefore given by 



(18) 



5r£ sym 



+ 



dn ~ % 
del 



( On 

- 2 9o(ph q B 



f {f+ + 5n+ + f+Sn') 



With the help of Eqs. (|T7|) and (fT9|) the average magnetization becomes, 



(19) 



s 



<M> 



|fe 2 ^[iV + (0) + iV-( 0)] 
l + [iV+(0) + iV-(0)]/ a 



asyra i 



(20) 



where we have suppressed the spin indices 
compared with the one presented in 



15 



'or /, 



s,asym(sym) 



The expression of (M) may be 



16[ to see the difference between the unpolarized 



and polarized matter. Likewise, the magnetic susceptibility is found to be 



Xm 



s g DN y [7V+(0) + iV-(0)] 



V 2 y 1 + [N+(0) + N-(0)}f° sym 



(21) 



where g D is the angular averaged gyro magnetic ratio 15 



3. 



With the help of Eqs.fjSJ) and ( 12~TT) we express the magnetic susceptibility in terms of LPs 



as, 



Xm 



Xp 



2tt 2 



rasym 
JO 



1 rsym 
3/1 



(22) 



15 



16| 



Here, xp — 9±j^lN c (p~ffi + +pJfi~)/(87i 2 ) is the Pauli susceptibility [15|, ll6[. For unpolarized 
matter £ = 0, implying p~f = pj, fi + = \i~ and iV + (0) = iV~(0). From Eq. (l2"2l we get the 
well known result for magnetic susceptibility 

NcPf/i 



Xm 



Xp 



7T Z 



rasym 

Jo 



esym 



3 fl J 



(23) 



A. Susceptibility with bare propagator 



We have already mentioned that the individual LPs are infrared divergent when evaluated 
with the bare gluon propagator. But the combination (fQ Sym — ^fi V j is always finite and 
turns out to be 



rasym x rsym 
Jo ~ g/l 



d(co S e)(i-co S e)(f+t + f-r 

d(cos0)(l + cos0)(/+7 + /^ 



h-h 



(24) 



Using Eqs.©, © and Eq.§H§ we have 



9 



h 

h 



9 2 f 1 

36lp+ 2 4 2 



2m z Q - p + f 2 + 
1 



2m q p^ 2 
f ~" 3(£+ + m q ) 



(25) 



x 



36ejey 3pJ ^ (m q + ej")(m ff + e / ) 



x < - 



-m, 



2pfP f p f e + f {m q p f - 3m q p+ - 3p + f e f ) + m 2 q (p+ 2 - 3p~fP f + p f ' 



,p+e f (p+ - 3pj 



(m q p f -PfP f -pj e f -m q e f [pj +p f \) 



,+2„-2 



+2 



Pf 2£ f £ V + ( m 9 -PfPf)( £ fPf' 2 + m g[pT +Pf 2 }) 



+2 



.+2 



In ' m - 



PjPt ~ e l e l 



mi 



PfPf 



£ f £ f 



) }(26) 



To determine \m for various £, we insert Eg. f)24p in Eg. ( 12 2 j) where Ji and /2 are given by 
Egs.()25]) and (EB)). 



B. Susceptibility with HDL corrected propagator 



In this section we consider the screening effects due to HDL corrected propagator of the 
j;auge field The scattering amplitude can be written as 15] 



[r°°(Ps, P's')D 00 + T'(Ps, P's')Di 



(27) 



In the coulomb gauge, we have D 00 = A; and = (5ij — qiqj/q 2 )A t , where q = p — p' . A 



and At denote the longitudinal and transverse gluon propagators given by 



20] 



1 

— r > A * 

q z + m D 



1 



ql - q 2 



(28) 



For spin parallel (s = s') and anti-parallel (s = —s') interaction, A; and A t have the following 
form : 



A t (s = s' 
A t (s = -s 
and A t (s = s' 



p=p<=p± 



p=p},p'=p f 



p=p<=p± 



2p^ 2 (l — cos#) + m 2 D 



AAs 



p=p J f,p'=v f 



P~f + Pf ~~ 2p~fPf (1 — cos I 

1 

~2pf{l -cosfl) ' 
1 

2{m 2 — £/ £ / + P~fPj cos @) 



m 



D 



(29) 



(30) 
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The matrix element given by Eq. (1271) can be calculated easily with OGE. We find that 



4m^ 



2p p' Q -P.p' + m 2 + (m 2 -P- P')(2a b - a ■ b) 



+2a p' (P ■ b) - 2p oP ' (a ■ b) + 2p b (a ■ P') - (P ■ b)(P' ■ a) (31) 



and 



T ij {Ps,P's') = Tr[fp(P,s)-f j p(P',s')} 

= {(1 - a ■ h )v l V ,j + (m 2 q - P ■ P')aW + (a • P')p^ 

+ {b ■ P)p^J + g ij [(m 2 q - P ■ P')(l - a ■ b) - (P ■ b){P' ■ a)]}, (32) 

with a symbol a % W = a l V + b % a? . Here, the spin vector and b^ are given by 



s + 



m q (e p + m q ) 
p'(p' ■ s') 



m q (e p > + m q 



p ■ s 
m q 
p' ■ s' 

m n 



(33) 
(34) 



To evaluate the spin symmetric and spin anti-symmetric combination of LPs, we need to 
calculate the scattering amplitudes both for spin non-flip (s = s') and spin flip (s = —s') 
interactions. The traces relevant for the longitudinal gluon exchange are given by, 



T 
In 



1 



oo 



Qm q (m q + £f) 2 



\2m\ + \2m\e) + Qm q e + f p + f 2 (l + cos 



+2, 



6m^ 2 (2 + cos 9) + p+ 4 (2 + 3 cos 



+4/ 



'00 



+2 -2 
Pf Pf 



Qm 2 q (m q + ef)(m q + e f ) 
Similarly, the coefficient of A t turns out to be 



(35) 
(36) 



%, x (<$*■ 



q % q> , 



v? 



6m q p^ 2 + 2p J j 2 e j 



Qm? q (m q + e + f 
+2m 2 q e+{4 + 3 cos 9) + m 3 q {8 + 3 cos 9) 



(37) 
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q l q J 



6m^(m q + e~f)(m g + £ f )(p~f 2 +p f 2 - 2p^p f cosi 



+m q (m q + 6f) - p f z {m q [p + f 2 +p/] +p/^ +P/^ £ /) 



-2_+ 



+2(p+ 2 - 3p+ Pf + p f 2 ){m q + e f ){m 2 q -ep] 



-PfPf 



cos 9 



2pf{m q \pf + p^] + pfe) + pjr 2 ^) 



,+2 c 



+m q (2m q [2p~j 2 - Sp^p f + 2j^ 



—2 +2 i +2 —2 
Pf £ f + PJ £ f 

H-fQ^+2 



+m q e f [hp j - QPfPf + Sp f ] + m g £^ [3p^ - QpjP f + 5p 
+3eJ"eJ [p + f 2 + Pj 2 - 2pjpJ cos 6}) 



(38) 



Using Eqs. fl271) - fl30l) and fl35l) - fl38l) one can easily calculate the required combination 
{fo Sym — jjfi V j to evaluate the magnetic susceptibility. Inserting this particular com- 



bination of f and f\ in Eq. f[2~2l we get Xm- To determine xm, we need to evaluate first /i + 
and n~ . This can be done by adopting the procedure outlined in Ref. 12]. With these, we 
can estimate xm numerically for the polarized and unpolarized matter at various densities. 
The corresponding results are discussed below. 



In Fig J2] we plot the magnetic susceptibility of cold and dense unpolarized quark matter as 
a function of Fermi momentum. It is observed that, upon inclusion of the screening effects, 
the divergence move towards lower densities. This is consistent with what one obtains for 
unpolled m atte r H Q. Such shifts toward .owe, density ar e expected as we hnow, that 



15 



16|). Moreover, this 



the screening effect weakens the Fock exchange interaction (See Ref. 
divergence is related to the magnetic phase transition of quark matter which shows up when 
the square bracketed term in Eq. (l22]) vanishes. As noted in [s| and also in [l2], this density 



approximately corresponds to the critical density for para-ferro phase transition. For the 
numerical estimation, we take a c = g 2 /4:ir = 2.2 and m q = 300MeV js], Q, Q, Q, 



In FigEJ the density dependence of magnetic susceptibility both for unpolarized and 
polarized matter has been shown. We see that the magnetic susceptibility diverges at some 
critical density which increases with increasing £. It is apparent from the figure that, if the 
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FIG. 2: Density dependence of magnetic susceptibility. Screening effects (solid line) are compare 
with the simple OGE case (dashed line) for unpolarized quark matter. 
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FIG. 3: Magnetic susceptibility vs Fermi momentum using screened gluon mass for unpolarized 
and complete polarized quark matter. 



critical density for para-ferro phase transition becomes lower than the critical density for 
the magnetic transition, the latter cannot take place. Thus, we conclude, that the magnetic 
transition depends on the critical density of para-ferro phase transition. 

In FigJU we show £ dependence of the magnetization for various densities. Note that 
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FIG. 4: Variation of magnetization with £ for a magnetic field 5 = 10 7 G. Solid, dashed and 
dash-dotted lines represent O.lfm -3 , 0.15fm -3 and 0.25fm -3 respectively. 



the divergences appear at higher £ for larger density. Here the magnetic dipole moments of 
the quarks are taken to be: fi u = 1.852/ijv, H>d = — 0.972/iiv and fi s = — 0.581/ijv, where the 
nuclear magneton /xjy = 3.152 x 10~ 14 MeV/ Tesla ^21 ] . 



IV. SUMMARY AND CONCLUSION 

In this work, we calculate dimensionless LPs F SJ ( m and i^ ym for dense quark matter. 
These are then used to calculate magnetic susceptibility and magnetization of degenerate 
quark matter and the results are found to be consistent with previous calculations in the 
appropriate limits. The qualitative behavior of the FLPs as a function of density is also 
found to be very similar to those of nuclear matter having isospin asymmetry. 

We observe that xm is free of all the infrared divergences even in the massless gluon limit. 
It is, however, numerically sensitive to the Debye mass. It is shown that the critical density 
for the magnetic transition in polarized matter is higher than that of the unpolarized one. 
The divergence and sign change of the magnetic susceptibility signal the magnetic instability 
of the ferromagnetic phase. 
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